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Charge-neutral excitons in semiconductor quantum dots have a small finite energy separation caused by 
the anisotropic exchange splitting. Coherent excitation of neutral excitons will generally excite both exciton 
components, unless the excitation is parallel to one of the dipole axes. We present a polaron master equation 
model to describe two-exciton pumping using a coherent continuous wave pump field in the presence of a 
realistic anisotropic exchange splitting. We predict a five-peak incoherent spectrum, thus generalizing the 
Mollow triplet to become a Mollow quintuplet. We experimentally confirm such spectral quintuplets for 
In(Ga)As quantum dots and obtain very good agreement with theory. © 2013 Optical Society of America 
OCIS codes: 130.5590, 270.0270, 300.0300. 



Quantum light-matter interaction in semiconductor 
nanostructures is a topic of considerable interest. An 
understanding of the fundamental light-matter inter- 
actions is important for future quantum devices, with 
applications ranging from single photon emitters [1, 2] 
entanglement-based sources [3,4]. Due to efficient semi- 
conductor growth technology and natural scalability, 
semiconductor quantum bits (QDs) offer an excellent 
opportunity to explore rich light-matter interaction in 
a solid state environment, and they can function as solid 
state quantum bits (qubits). In the presence of high- field 
coherent optical pumping, it is well known that the flu- 
orescence spectrum of a driven two-level system devel- 
ops a symmetric "Mollow triplet" structure [5], where 
the two outer sidebands are split by the drive's Rabi 
frequency. While driven QDs, also referred to as "ar- 
tifial atoms," do show behaviour that is similar to a 
driven two- level atom [6,7], due to the coupling of the 
exciton to a phonon reservoir [8-13], both the position 
and broadening of the Mollow triplets can change signif- 
icantly [14, 15]. It is also well known that real QDs have 
many exciton states, and, e.g., charge- neutral excitons 
are split by a small anisotropic exchange energy. 

In this Letter we develop a polaron master equation to 
include two neutral excitons and model the ensuing inco- 
herent spectrum. We generalize the well known Mollow 
triplet to a Mollow quintuplet regime, which is caused 
by a sum of two separate Mollow triplets with a similar 
central resonance. A useful analytical expression for the 
incoherent spectrum is given and compared to experi- 
ments, where we find very good qualitative agreement. 
Figure 1 shows a schematic of the excitation geometry. 

Neglecting QD zero-phonon-line (ZPL) de- 
cay mechanisms, and working in a rotating 
frame with respect to the pump laser frequency 
cjl, the two-exciton Hamiltonian is given by 
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(j = 1, 2) are the laser-exciton detunings, 6jq are the 
annihilation operators of the phonon reservoirs, a^, a~ 
are the Pauli operators of the jth exciton, Ajq represent 
the exciton-phonon coupling strength, and are the 
Rabi frequencies of the continuous wave (cw) pump. 
Adopting a polaron transformation R' = e^He~^ , 

with S = E^q^/^7^(^]q - ^jq), thc CXCitOU- 

phonon coupling can be taken into account non- 
perturbatively [11, 12]. To second-order in the 
polaron-transformed exciton-phonon interaction, 
the polaron master equation can be written as 

^ dT{[Hi,Hi{-T)p{t)] + H.c), where the upper 
time limit of the time integral, t oo, since the 
relaxation time of acoustic phonon bath is very fast 
(a few ps). The polaron- modified system Hamiltonian, 

where the renormalized Rabi frequency is given by Q^j = 
^j{Bj), with {Bj 
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(a) Planar Sample Structure 




(b) QD Side-Excitation 



^QD 



Fig. 1. Schematic of the QD under continuous wave (cw) 
pumping with a field that is polarized along, e.g., = 
^ in a weakly coupled planar cavity. This QD can be 
excited with electric dipole along the x or y axis. 
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which depends on the temperature of the phonon 
bath. We chose a phonon spectral function that 
accounts for exciton-LA interactions from the defor- 
mation potential, Jj{uj) = a^^cj^ exp(— a;^/2a;^.) [16], 
where 00^. is the phonon cutoff frequency and ap. 
characterizes the strength of the exciton-phonon in- 
teraction. We also include Lindblad scattering terms, 
C[D] = {DpD^ - D^Dp) + H.c, which describe 
ZPL radiative decay and ZPL pure dephasing [17]. 
For convenience, we will include the polaron shift 
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For Rabi fields much less than uji). , one can derive an 
effective phonon master equation [13,14], 



j 



E 



_ph 
2 



Ei^^[-7]-Erpl,>5, (1) 



where S = {gJ paj + a~pa~). This model is ap- 
plicable to two separate QDs or two excitons from 
the same dot, and we will focus on the later here to 
model realistic single QDs with driven neutral excitons. 
The phonon-induced scattering rates are derived as fol- 



lows [13,14]: Fpf^ = ^Re[f^ dre^'^^^ {e' 
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the one exciton limit, the above model has been used to 
describe the spectra for driven In(Ga)As QDs [14], where 
phonon-mediated processes are also found to play a sig- 
nificant role. Consistent with these experiments, we will 
use ujbi = ^62 = 1 meV and ap^ = ap^ = 0.15/(27r)^ ps^. 

The emission from the QD is experimentally detected 
through a planar cavity (which spatially separates the 
pump field from the emitted spectrum), and the inco- 
herent spectrum of the system is given by [14] S{r^uj) = 
F{r)S{uj). Here F{r) is a geometrical factor, and S{uj) 
can be calculated by the quantum regression theorem. 
From our two-exciton polaron master equation, the an- 
alytical spectrum is calculated to be 
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^'y[2{i5u-%,^% and E^{u:) = 
The steady-state functions are 

(Sapa-),, = i [l + (a|)ss-2(a-)ss(a+)ss], 
{SajSaiU = and hj{0) ^ 



{ScrjScr^)ss = -(crj)ss [1 + (cr|)ss], where the steady- 
state inversion and polarization are given by (cr|)ss = 



" (7pop 



2r;o/bpop- 



"'j(7cd+7poi)-| 



and (a- 



0.6 



0.4 







(a) 


1 9 = t:/4 




{ A; 




1 e = TT/8 


. . rtA 


\ A A - 






AA J 






-50 50 -50 50 

u; — ujl (jiieV) uo — uol (a^gV) 

Fig. 2. (Color online) Incoherent spectrum S(uj) with 
(blue, solid) and without (red, dashed) phonon scat- 
tering. The Rabi frequencies are Vt'i — ft' cos 6 and 
Q2 = ^^'sin6>, and 71,2 = 7i 2 = ^ /xeV, S12 = — 12 /ieV, 
and ai,2 = 1- (a) ft' = 30 /ieV, with different 0. (b) 
6> = 7r/8, with different Q' . 
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lation decay rates are defined through ^j^^i 



7i + 7j + + T^i and 7^„p = 7i + + ■ Fi- 
nally, aj are phenomenal scaling terms that adjust the 
strength of the exciton emission through cavity decay. 

Due to Coulomb interactions and orientational in- 
homogeneities of the QDs, there are two neutral ex- 
citons polarized along the x and y axes, respectively 
[c.f. Fig. 1]. Thus for various pump polarization angles, 
^, there will be different Rabi frequencies, Q[ = Q[ cosO 
and Q2 — ^5 ^^o excitons. Here Vt[ = {B) 0!^ 

is the maximum Rabi frequency. For clarify, we will first 
assume that Q.'^ = = (B) = Q\ though in gen- 
eral these will differ because the dipole moments will 
in general be different. Figure 2(a) shows an example 
calculation of the incoherent spectrum for various ^, 
using Ai = — A2 = —6 peY. We assume a phonon 
bath temperature of T = 6 K, and use parameters for 
In(Ga)As QDs [18]. The main parameters are given in 
figure caption. The red dashed lines shown in the same 
figure are the results obtained without phonon scatte- 
ring (apart from ZPL decay), which show that phonon 
scattering plays a qualitatively important role; in partic- 
ular, we see phonon-mediated spectral broadening and 
phonon renormalization induced by phonon scattering. 
For ^ = 0, only one exciton is excited, and the spec- 
trum shows the expected Mollow triplet; note that the 
spectrum is asymmetric since the pump field is off reso- 
nance. As increases, then the other orthogonal exciton 
is gradually excited, as can be clearly seen at 6> = tt/S — 
where the Mollow triplet evolves into a Mollow quintu- 
plet. At = 7r/4, the quintuplets merge into a triplet 
again if the two excitons share symmetric parameters. 
Figure 2(b) shows the spectrum at ^ = | for various 
Rabi fields (I^q = 20, 30, 40 peV from bottom to top). 
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Fig. 3. (Color online) Experimental data for a driven 
In(Ga) As QD and corresponding theoretical simulations, 
both as a phonon bath temperature of 6K. (a) Experi- 
mental results with pump powers 50 fiW (dark, lower) 
and 100 /iW (red, upper), (b) Theoretical fits, with the 
following parameters: 71 = 0.55 /ieV, 72 = 0.85 /ieV, 
j[ = 5.93 /ieV, 7^ = 7.13 /ieV, S12 = 10 fieV. The Rabi 
frequencies are Q[ = 13.66/>/2 /ieV, ^2 — 21.96/>/2 /ieV 
(dark, lower), and Q.'^ = 13.66 /ieV, = 21.96 /ieV (red, 
upper); we also choose ai = 1 and 0^2 = 0.4. 



which show how the spectral quintuplets evolves with 
pump strength. From the analytical spectrum, it can be 
found that there will be six peaks for the incoherent spec- 
trum in general. However the central peak (for both ex- 
citons) will always lie at the laser pump frequency and 
show one peak, so there are actually five peaks. 

Next we turn our attention to experiments. The planar 
sample under investigation is grown by metal-organic va- 
por epitaxy. Self- assembled In(Ga) As QDs are embedded 
in a GaAs A-cavity, sandwiched between 29 (4) periods of 
A/4-thick AlAs/GaAs layers as the bottom (top) distrib- 
uted Bragg refiectors. For the investigations, the sample 
is kept in a Helium fiow cryostat providing high tempera- 
ture stability T = 6 ± 0.5 K. Laser stray- light suppres- 
sion is achieved by use of an orthogonal geometry be- 
tween QD excitation and detection in combination with 
polarization suppression and spatial filtering via a pin- 
hole. Resonant QD excitation is achieved by a narrow- 
band (~ 500 kHz) continuous- wave (cw) Ti: Sapphire 
ring laser. For high-resolution spectroscopy (HRPL) of 
micro-photoluminescence (/i-PL) we employ a scanning 
Fabry-Perot interferometer with A£^™^^ < 1 jieV as 
described earlier [14,18]. 

Figures 3(a)-(b) show the experimental results and 
theoretical simulations, respectively, for two different 
pump fields. Due to the fact the experimental central 
peak includes contributions from both coherent and in- 
coherent, we focus on reproducing the experimental data 
around the Mollow sidebands. Figure 3(a) is the ex- 
perimental results under pump strength 50 jiSSf (dark, 
lower) and 100 /iW (red, upper). The dark (inner) ar- 
rows and red (outer) arrows show the inner sideband 



and outer sideband of the quintuplet, respectively. Fig- 
ure 3(b) shows the theoretical simulations with fitting 
parameters given in the caption, which are consistent 
with earlier experiments [14, 18]. We obtain a very good 
qualitative agreement. Moreover, we have observed these 
spectral quintuplets from many of our In(Ga)As QDs. 

In summary, we have introduced a theory to describe 
resonance fiuorescence of coherently excited QD neutral 
excitons and predict the emergence of a spectral quin- 
tuplet. The theory uses a polaron master equation from 
which an analytical spectrum is derived and presented. 
Using data for In(Ga)As QDs, we obtain excellent agree- 
ment between theory and experiment. 

This work was supported by the National Sciences and 
Engineering Research Council of Canada, the Carl-Zeiss- 
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